Anomalous structure and dynamics of the Gaussian-core fluid 
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It is known that there are thermodynamic states for which the Gaussian-core fluid displays anoma- 
lous properties such as expansion upon isobaric cooling (density anomaly) and increased single- 
particle mobility upon isothermal compression (self-difltusivity anomaly). Here, we investigate how 
temperature and density afltect its short-range translational structural order, as characterized by the 
two-body excess entropy. We find that there is a wide range of conditions for which the short-range 
translational order of the Gaussian-core fluid decreases upon isothermal compression (structural 
order anomaly). As we show, the origin of the structural anomaly is qualitatively similar to that 
of other anomalous fluids (e.g., water or colloids with short-range attractions) and is connected to 
how compression affects static correlations at different length scales. Interestingly, we flnd that 
the self-diffusivity of the Gaussian-core fluid obeys a scaling relationship with the two-body excess 
entropy that is very similar to the one observed for a variety of simple liquids. One consequence of 
this relationship is that the state points for which structural, self-diffusivity, and density anomalies 
of the Gaussian-core fluid occur appear as cascading regions on the temperature-density plane, a 
phenomenon observed earlier for models of waterlike fluids. There are, however, key differences 
between the anomalies of Gaussian-core and waterlike ffuids, and we discuss how those can be qual- 
itatively understood by considering the respective interparticle potentials of these models. Finally, 
we note that the self-diffusivity of the Gaussian-core fluid obeys different scaling laws depending 
on whether the two-body or total excess entropy is considered. This finding, which deserves more 
comprehensive future study, appears to underscore the significance of higher-body correlations for 
the behavior of fluids with bounded interactions. 



Model pair potentials that describe the interactions be- 
tween particles of simple atomic and molecular liquids 
generally have steeply repulsive, short-range components 
that diverge as two particles approach one another, qual- 
itatively capturing the increase in energy that accom- 
panies overlap of their electronic clouds. However, the 
effective pair interactions between the centers of mass of 
two polymer chains in solution (or molecular aggregates 
such as micelles) are often more accurately represented by 
softer, bounded potentials [ij, since they describe "par- 
ticles" that are inherently interpenetrable. A simple ex- 
ample of a bounded interaction is given by the Gaussian- 
core (GC) model 0. It is defined by a pair potential of 
the following form, (j)(r) = eexp[— (r/cr)^], where r is the 
distance between particle centers, and e and a charac- 
terize the amplitude and width of the interaction profile, 
respectively. Understanding the differences between the 
collective properties of penetrable particles with bounded 
effective potentials and those with steeply repulsive inter- 
actions is a basic problem in the study of soft condensed 
matter (see, e.g., H, 

The behavior of the GC fluid is indeed unusual when 
compared to most atomic or molecular fluids. For cxam- 
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pie, there are state points for which its thermal expan- 
sion coefficient is negative (density anomaly) [6i] and for 
which its single-particle dynamics increase upon isother- 
mal compression (diffusivity anomaly) 0]. At suffi- 
ciently low temperature, the fluid also exhibits re-entrant 
melting behavior [1, H, Anomalous changes of self- 
diffusivity with density or temperature are known to oc- 
cur in other model systems that qualitatively mimic ei- 
ther the behaviors of (i) liquid water or (ii) suspensions 
of colloids with short-range interparticle attractions, re- 
spectively. A series of theoretical studies have shown 
that those changes for both systems of type (i) [1, ITol . 
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type (ii) ^Q, |23, [231 are strongly correlated to anomalous 
trends in the short-range translational structural order 
of the respective fluids. Moreover, an analysis of re- 
cent Brownian dynamics simulations of the GC fluid [l^ 
suggested a qualitative (but not quantitative) connection 
between the anomalous density dependencies of the self- 
diffusivity and the pair correlation function. Here, we 
carry out an extensive molecular simulation study of the 
GC fluid to investigate whether a quantitative link be- 
tween its dynamics and its short-range translational or- 
der can be established. 

To compute the properties of the GC model, we per- 
form molecular dynamics simulations in the microconan- 
ical ensemble using N ~ 1000 particles in a cubic cell 
with periodic boundary conditions. For simplicity of 
notation, we implicitly non-dimensionalize all quantities 
by appropriate combinations of the characteristic length 
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scale Ic = (T and time scale Tc = ^/ma^Je, where m is 
the mass of a particle. We integrate the equations of mo- 
tion using the velocity- Verlet method [l^l with time step 
St = 0.05, and we cut the GC potential at Tcut = 3.2. We 
calculate sclf-diffusivity by fitting the long-time (t 3> 1) 
behavior of the mean-squared displacement {Sr'^{t)) to 
the Einstein formula 6Dt = {6r'^{t)). We investigate tem- 
peratures between T = 0.01 and 1.5 and number densities 
between p = 0.01 and 2.0 [1^, parameters that span the 
part of phase space where the aforementioned anomalous 
behavior is known to occur for the GC fluid Q. 

Figure [Tal displays the self-diffusivity D of the GC fluid 
versus density p for several isotherms. Like in simpler 
fluids, compressing the low density GC fluid initially de- 
creases D. However, as noted previously f?\, compression 
can eventually lead to an anomalous increase in D at suf- 
ficiently high particle density. In fact, for a given T, one 
can view the density at which the minimum in D occurs 
as a boundary between regions of "normal" and anoma- 
lous dynamic behavior. This boundar y is approximately 
indicated by the green line in Fig. [Ta| [29| . 

To address the question of whether this anomalous dy- 
namic behavior is related to structural order, we first 
investigate the behavior of the two-body contribution to 
the excess entropy (relative to an ideal gas at the same 
T and p), defined as 

nOC 

S2 = -2ttp / r^g{r)\ng{r) - [g{r) - l]}dr, (1) 

where g{r) is the pair correlation function (PCF). The 
quantity —S2 effectively characterizes the degree of pair 
translational order present in the fluid [32| . FigurefTbldis- 
plays — S2 versus p for several isotherms. Like in simple 
fluids, compressing the low density GC fluid initially in- 
creases structural order. However, further increasing the 
density eventually leads to a qualitative change in behav- 
ior. Specifically, the GC fluid begins to lose structural or- 
der upon compression at sufficiently high density, a trend 
that docs not generally occur in atomic fluids [Sll, [s^ . 
Thus, for a particular T, the density corresponding to the 
maximum in —S2 can similarly be viewed as a boundary 
between state points of "normal" and anomalous struc- 
tural order. This boundary is approximately indicated 
by the green line in Fig. [lb] [1^ . 

Together, Figures [la] and Jb] suggest a strong correla- 
tion between dynamics and structural order for the GC 
fluid. Given the recent results from computer simulations 
[11, [11, 111, [Ipjpjp [11, [11, 111 [11, [11, IS and 
experiments fij, [4j which show a clear connection be- 
tween dynamics and excess entropy for a wide variety of 
systems, this correlation may not be particularly surpris- 
ing. In fact, an approximate scaling law relating trans- 
port coefficients and excess entropy per particle, s°^, was 
first noted by Rosenfeld [H, [s^ . Specifically, Rosenfeld 
observed that the following scaled form of self-diffusivity, 
Dft. = Dp^/^T~^/'^ ^ is a quasi-universal function of s^'^ 
for a variety of models of simple atomic fluids. In other 
words, the dependencies of D on T and p can essentially 



be reduced to how these variables affect s'^^. For simple 
liquid systems at low to intermediate densities, most of 
the excess entropy {s°^) comes from the two-body con- 
tribution (S2) 31], and thus the quasi-universal excess 
entropy scaling indicates that their transport coefficients 
are closely linked to g{r) (see also, e.g., p5j). 

For atomic fluids, the scaling law relating D and s"'^ 
closely tracks different functional forms as —s°^ goes from 
lower to higher values. For very low — (i.e., a dilute 
gas), the following power-law form follows from Enskog 
theory for D and a second-virial approximation for 
(see |35|), 



l/3rTn-l/2 



Dp^'-'T 



a(-s°") 



,ox\/3 



(2a) 



where a, (3 are constants. At higher — s°'^ (e.g., a fluid 
near the freezing transition), the following exponential 
relation has been empirically observed [351] 
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where A and B are constants. 

Figure [Tc] displays the scaled self-diffusion coefficient 
Dji versus the two-body estimate of the excess entropy 
—S2- Note that the dependencies of the Rosenfeld self- 
diffusivity on T and p for the GC fluid also collapse 
onto a single master curve when plotted versus —S2- For 
comparison, Eq. ((^J) with parameters that fit the sim- 
ulation data for a variety of soft-sphere model systems 
{a = 0.37, /3 = -2/3 [ll and A = 0.58, B = 0.78 [H) 
is also shown in Fig. [Tc] The main point to note here 
is that, despite small quantitative differences, the form 
of the relationship between the diffusivity and excess en- 
tropy is virtually the same, whether one considers the 
anomalous GC fluid or a model for a simple atomic fluid. 
In other words, despite the anomalous dependency of D 
upon p for the GC fluid, the relationship between D and 
S2 is perfectly normal. 

Given this basic link between anomalous dynamics and 
structure, we now investigate the molecular origins of the 
trends in —S2- To this end, we study the PCF and the 
related cumulative order integral 1^^ (r) defined as [2^ 



/..(r) = 27rp 



r'^{g{r')\ng{r') - [g{r') - l]}dr' . (3) 



Note that one recovers —S2 from this integral in the large 
r limit. The quantity Is^i''') characterizes the average 
amount of translational ordering on length scales smaller 
than r surrounding a particle. Figs. Ildl - |g] show the 
behavior oi g{r) and Ig^ [r) for selected state points along 
the low-temperature T = 0.025 isotherm (red curve in 
Fig. Ilb[) . They divide the behavior into two qualitatively 
different regions: (1) the normal increase of —S2 with p 
at low densities [p < 0.2, Figs, [id] and |ej and (2) the 
anomalous decrease of —S2 with p at high densities [p > 
0.2, Figs. [13 and [g. 

Inspection of Fig. [id] shows that compression of the 
low-density GC fluid leads to an increased degree of 
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FIG. 1: (Color online) (a) Self-diffusivity D and (b) structural order metric — S2 versus density p at temperatures 
T = 0.01,0.013,0.018,0.025,0.035,0.05,0.07,0.1,0.15,0.2,0.3,0.4,0.5,0.7,1.0,1.2, and 1.5 for the GC fluid. Arrows indicate 
increasing temperature. Solid green lines in (a) and (b) correspond to minima in D and maximum in — S2 , respectively, 
(c) Rosenfeld scaled difTusivity versus — S2 at all densities and temperatures studied. The solid line represents the empiri- 
cal Rosenfeld scaling law 0.58exp(0.78s2) and the dashed curve represents the low density behavior expected for soft spheres 
0.37(— S2)~^''^. (Lower right panel) Pair correlation function g{r) and cumulative order integral Is2{i') along the isotherm 
T — 0.025 [red curve in (b)] for two density ranges [(d) and (e)] p < 0.2 [below the maximum in — S2(p)] and [(f) and (g)] 
p > 0.2 [above the maximum in — .S2(p)]. In (d-g), arrows indicate increasing p. 



translational order, reflected by the formation of distinct 
first, second, and more distant coordination shells in g{r). 
This increase in "packing" order around the particles (see 
Fig. [Te)) is also observed in hard-sphere fluids and simple 
atomic systems, and it is a reflection of the fact that local 
interparticle correlations necessarily build up as density 
increases in order for particles to avoid overlap with one 
another. At sufficiently low pressures (low densities) and 
low temperature, the effective repulsive "core" of the GC 
potential acts like a hard-sphere interaction, and that is 
the primary reason that the structural trends of the GC 
fiuid under these conditions mimic those of simple atomic 
systems. 

Further compression of the GC fluid [region (2), p > 



0.2] leads to an anomalous decrease in structural order 
[Fig. lib] . In this higher pressure, higher density region of 
phase space, the repulsion of the GC interparticle poten- 
tial is too soft to effectively prevent interparticle overlap. 
In fact, because the system only pays a finite potential 
energy penalty for each overlap, and since avoiding over- 
laps carries a substantial entropic penalty, the system 
evolves toward a more uniform average structure at suffi- 
ciently high density. In the limit of infinitely high density, 
the particle centers of the GC fluid will adopt a Poisson 
distribution (i.e., a high-density ideal gas) As can 
be seen in Figure [T3 these compression induced changes 
manifest as a shifting in, flattening, and broadening of 
the coordination shells in g{r). Figure [Tg| shows that 
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FIG. 2; (Color online) (a) Boundaries of the structural, diffusivity, and density anomalies in the temperature-density plane, as 
calculated from simulation data. The inset shows the boundary of the region of density anomalies on an expanded scale. The 
symbols are the locations of the extrema, and arrows are in the direction of the anomalous regions, (b) Schematic representation 
of anomalous regions. 



while finite penetrability gives rise to a small increase in 
order at small r (inside the core) , the overall behavior is 
dominated by the corresponding disordering that occurs 
at larger r. 

The microscopic origins of the structural behavior of 
the GC fluid share a qualitative similiarity with those 
of the other anomalous fluids discussed earlier (liquid 
water and suspensions of colloids with short-range at- 
tractions) [2^. In the case of models that capture wa- 
terlike behavior, compression of the low temperature 
fluid can result in an anomalous translational disorder- 
ing 0, M, nil, M, [H, [H, H m, [13, m, Si, 113, S, m. 

This decrease in structural order is primarily due to a 
shifting in, flattenin g, a nd broadening of the second co- 
ordination shell (20l. |23||). which is accompanied by the 
penetration of a fifth water molecule into the periphery 
of the first shell of otherwise tetrahedrally-coordinated 
molecules (see, e.g., ^3)- ^^e case of colloids with 
short-range attractions, increasing the strength of the 
attractive interactions (relative to fceT') of the concen- 
trated fluid can also give rise to an anomalous disorder- 
ing [l^ [2^ . This disordering is due to the fact that the 
short-range attractions favor near contact conflgurations 
of the particles. This disrupts the structure of the second 
and more distant coordination shells that otherwise natu- 
rally form when structure is determined primarily by the 
repulsive interactions (i.e., as in hard-sphere or simple 
atomic fluids) 23J. Thus, in each of the above systems, 
there are two effects that are linked: (i) an increase in 
correlations at short length scales and (ii) an associated 
decrease in structural order at larger length scales. The 
latter effect dominates in these fluids, thus giving rise to 
their structurally anomalous behavior. 



It is also interesting to ask about the relative placement 
of the boundaries for structural, diffusivity, and density 
anomalies on the temperature-density plane. Fig. I2al 
shows those boundaries for the GC fluid as determined 
from our simulation data, and Fig. [2b] provides a ideal- 
ized schematic based on that data. The most striking 
feature is that the state points which exhibit the density 
anomaly (negative thermal expansion coefficient) are a 
subset of those associated with the diffusivity anomaly, 
which, in turn, are a subset of those associated with the 
structural anomaly. This "cascade of anomalies" is qual- 
itatively similar to that originally noticed for a molec- 
ular model of water by Errington and Debenedetti 
and more recently observed in a wide class of model sys- 
tems that exhibit waterlike behavior [lol. [Til , [l^ . [Tsl [3, 
m, [13, m, m, M, mi, m, 113, m, m, 131. one qualita- 
tive difference between the GC fluid and water is that, 
while water returns to normal structural and diffusiv- 
ity trends at sufficiently high density, the GC fluid does 
not. This difference follows from the fact that water has 
steeply repulsive short-range interactions, which domi- 
nate its physics at sufficiently high densities (pressures). 
The repulsive core is, of course, absent in bounded po- 
tentials like the GC model, and thus those systems never 
"return" to hard-sphere-like structural and dynamic be- 
havior at high density. However, like water, the GC fluid 
does return to normal (i.e., positive) thermal expansion 
behavior at sufficiently high density. This is consistent 
with the fact that the GC fluid's equation of state ap- 
proaches a mean-field description in that limit Fi- 
nally, we note that de Oliveira et al. (s^l have created a 
simple model for water with a spherically-symmetric pair 
potential represented by a sum of a Lennard- Jones contri- 
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bution (which contains the necessary repulsive core) and 
a Gaussian-core interaction. Given the discussion above, 
it should not be surprising that this generic type of hy- 
brid model can qualitatively reproduce water's density, 
diffusivity, and structural anomalies i4^, \^ . 

It is also interesting to note that the state points where 
models of colloidal fluids with short-range attractions 
show diffusivity anomalies are similarly a subset of the 
state points where they exhibit structural anomalies [1^ . 
As has now been discussed extensively for models of wa- 
ter [l7|Land for models of colloids with short-range attrac- 
tions 25], the "cascade of anomalies" has been shown to 
be qualitatively consistent with the behavior of the ex- 
cess entropy, the empirical relationship between diffusiv- 
ity and excess entropy, and the rigorous thermodynamic 
link between the excess entropy and the thermal expan- 
sion coefficient. 

Thus far, we have focused exclusively on the two-body 
contribution to the excess entropy (52) of the GC fluid 
rather than the excess entropy itself (s®''). The ratio- 
nale for doing so is simple. The two-body quantity is 
a particularly convenient structural measure to compute 
and analyze, given its straightforward connection to g{r). 
Furthermore, as discussed above, S2 is known to closely 
approximate for simple fluids at low to intermediate 
density [3l|. However, the GC fluid has structural be- 
havior that is very different from that of simple atomic 
fluids. In fact, it is known that there are significant 
differences between S2 and for this system, due to 
non-canceling effects from higher-order static correlations 
(3-body, 4-body, etc.) 51]. As an example, along the 
T = 0.025 isotherm highlighted above, the ratio s°'^/s2 
evaluates to approximately 1.7 at p = 0.3 (where —s°^ 
reaches a maximum) and 3.0 at p — 1.0. Thus, as a 
final point, we present in Fig. [3a| the relationship be- 
tween the Rosenfeld diffusivity, Dj^ — Dp^^^T^^/'^ , and 
—5°'' for all state points considered in this study of the 
GC fluid. We calculate in this study using grand- 
canonical transition-matrix Monte Carlo (GC-TMMC) 
simulations [13, [H, l53|, IsJ] . 

Note that, for a given density, the relationship be- 
tween Dn and —s°^ is qualitatively similar to that of 
atomic fluids and hence also to the relationship between 
Dr and — S2 for the GC fluid shown in Fig. [Tel How- 
ever, the Dfi data at different densities do not collapse 
onto a single curve when plotted versus —s'^'^ (contrast 
with Fig. [Tc)) . As expected, the low-density data are 
indeed accurately represented by Eq. ^ with parame- 
ters that describe soft-sphere model systems. However, 
as density is increased, Dr becomes increasingly under- 
estimated by the soft-sphere relation. Interestingly, as 
is shown in Fig. I3b[ the quantity DT~^/'^ does approxi- 
mately collapse for all densities when plotted versus — s"'^ 
for the GC fluid. At present, it is not known why this 
alternative scaling holds for this system or whether it 
can be expected to describe the behaviors of other flu- 
ids with bounded potentials (e.g., the penetrable sphere 
model 3]). We are currently carrying out a systematic 
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FIG. 3: (a) Rosenfeld scaled diffusivity versus excess entropy 
—s"^^. The solid and dashed curves, as well as the symbols, 
are the same as those in Fig. [Tc] The arrow is in the direc- 
tion of increasing density, (b) Temperature scaled diffusivity 
j)y-i/2 Yerses — s™. The data represent all densities and 
temperatures considered in this study. 



investigation of the thermodynamic and dynamic behav- 
ior of a variety of bounded potential fluids to explore this 
issue, and we will report our findings in a future study. 

In conclusion, we present molecular simulation re- 
sults which demonstrate that the GC fluid exhibits cas- 
cading regions of density, self-diffusivity, and structural 
(two-body excess entropy) anomalies in the temperature- 
density plane. Interestingly, an appropriately normalized 
self-diffusivity of the GC fluid follows a scaling with the 
two-body excess entropy that is similar to that observed 
for simple atomic fluids, which do not exhibit anomalies. 
The significance of this relation is that the self-diffusivity 
anomaly can be viewed as simply a reflection of the struc- 
tural anomalies of the GC fluid. An analysis of the pair 
correlation function shows that the molecular origins of 
the structural anomalies of the GC fluid are easy to un- 
derstand and are qualitatively similar to those of liquid 
water and colloids with short-range attractions. Specifi- 
cally, changes in thermodynamic parameters give rise to 
an increase in pair correlations at short length scales and 
an associated disordering at larger length scales. For all 
of these fluids, the latter effect dominates, giving rise to 
the anomalous structural trends. Finally, we show that 
while the Rosenfeld diffusivity [Dr = Dp^l'^T'^l'^) is 
approximately a function of S2 alone for the GC fluid, a 
different combination {DT^^^^) collapses as a function of 
s'^^. Whether this latter observation will generally hold 
for other fluids with bounded potentials is currently an 
open question which deserves more careful study. 
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